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EVALUATION OF BOUNDARY- LAYER AND 
WAKE-SURVEY DATA- REDUCTION TECHNIQUES 
IN COMPRESSIBLE FLOWS 

By K. R. Czarnecki and Russell B. Sorrells HI 
Langley Research Center 

SUMMARY 

An analytical and experimental investigation has been made of boundary- layer and 
wake-survey data- reduction techniques in compressible flows where the local Mach num- 
ber just outside the boundary layer or wake at the measuring station does not match the 
average free- stream or normalizing value. The analytical study was limited to turbulent 
boundary layers with no heat transfer and to a Mach number range from 0 to 10. The 
experimental portion of the investigation was confined to zero-heat-transfer, turbulent- 
wake surveys made behind a swept wing at free-stream Mach numbers of 1.61 and 2.01. 

Results of the analytical boundary-layer calculations indicate that if simplified data- 
reduction techniques that neglect the deviation of the local Mach number from the average 
free-stream or normalizing value are utilized, sizable differences in momentum, dis- 
placement, and velocity thicknesses are incurred over much of the Mach number range for 
relatively small deviations in the local-flow characteristics or parameters. In general, 
use of the free-stream Mach number as a base both for determining velocity defects due 
to skin friction and for normalizing the thickness parameters, will introduce larger dif- 
ferences than use of local Mach number as the base for both items and interpretation of 
results as being obtained under stream conditions. Experimental wake data indicate that 
because of increased wake widths and reduced viscous heating or compressibility effects 
in comparison with boundary layers of equal momentum deficiency, the use of the free- 
stream-flow technique generally results in proportionately larger differences in wake- 
data reduction. With proper reservations for basic differences between wakes and 
boundary layers, the trends in the differences of the experimental wake data appear to be 
in good agreement with the analytical calculations for boundary-layer results which are 
reduced by the same method. 


INTRODUCTION 

During the experimental investigation by means of wake surveys of three- 
dimensional -wing skin-friction drag at supersonic speeds (ref. 1) an approximate 



evaluation was made of the probable magnitude of the differences that might be incurred 
by the use of various methods of data reduction for the case where the flow conditions at 
the measuring station differ slightly from the desired reference conditions. This limited 
evaluation indicated the possibility of large differences under conditions where a cursory 
check of the accuracy of the basic wake- survey data would appear to justify the use of 
simplified data- reduction techniques. It thus appeared desirable to investigate this pos- 
sibility of large differences in calculated boundary- layer parameters in greater detail. 

For convenience, a basically analytical approach was chosen, which was, however, 
supplemented by a study of some experimental data. The analytical approach was limited 
to boundary layers inasmuch as the basic problem is pertinent in the reduction of 
boundary- layer data, and the theory for the development of wakes in compressible flows 
contains too many gaps (ref. 2) that cannot be handled analytically with any confidence. 
For further simplicity, the analysis was made for compressible-flow turbulent boundary 
layers with no heat transfer and a Prandtl number of 1, for boundary layers having 
power-law-type velocity profiles, and for isentropic free-stream flow. It should be 
recognized, nevertheless, that since similarities in the essential features of both laminar 
and turbulent boundary- layer flows exist, the trends derived from the turbulent boundary 
layers will also apply approximately to laminar boundary layers. Two types of approxi- 
mate data-reduction techniques were studied. The boundary- layer parameters analyzed 
were the momentum, displacement, and velocity thicknesses. The Mach number range 
for the calculations was from 0 to 10. In the experimental study, some of the wake- 
survey data of reference 1, taken at varying distances downstream of the trailing edge of 
a sweptback wing at Mach numbers of 1.61 and 2.01, were reduced to momentum, dis- 
placement, and velocity thicknesses by the same methods; the results were assessed in 
terms of relative scatter and deviation from the results computed by the more correct 
reference method, which is described subsequently. 

SYMBOLS 

a to g constants (defined by eqs. (33)) 

Cj section skin-friction coefficient (defined by eq. (38)) 

c s section of wing chord 

E difference function (defined by eqs. (13) to (18)) 

M Mach number 
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Mach number deviation ratio 
summation indices (see eqs. (33)) 
static pressure (appendix C only) 

velocity-profile power- law parameter (defined by eq. (25)) 
velocity in stream direction 

velocity transformation variable (defined by eq. (B7)) 
longitudinal distance downstream of wing trailing edge 
distance normal to surface 

ratio of specific heat at constant pressure to specific heat at constant volume 

total boundary-layer thickness 

boundary- layer velocity thickness, f (l - — W 

J 0 \ u r/ 


boundary-layer displacement thickness, 
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r\ 0 

boundary- layer momentum thickness, \ — (1 - v~^dy 

J 0 p co u oo\ u r/ 

mass density of air 


Subscripts: 
r reference 

6 local conditions outside boundary layer or wake except as noted in 

appendix C 

Su boundary- layer velocity thickness 

6* boundary- layer displacement thickness 
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boundary-layer momentum thickness 


°° free- stream conditions 

Superscripts: 

r reference method 

A local-flow method 

B free -stream flow method 


THEORY 

The basic problem, which is the subject of this paper, arises because the local 
stream conditions just outside a boundary layer at the measurement station are often not 
identical with the average or reference free- stream values owing to the existence of 
pressure gradients or other free- stream disturbances, and yet the differences are so 
small that there is a tendency to ignore them. For such conditions, it is assumed that the 
expressions for the losses within the boundary layer of momentum, mass flow, and 
velocity, respectively, can be determined from the local conditions as follows (see the 
shaded portions of sketches in fig. 1): 


r 6 

J 0 pu ( u « • ") dy 

(1) 

x W>« ■ pu ) dy 

(2) 

( u 5 - u ) dy 

(3) 


where ug is the local velocity just outside the boundary layer. Appendix A justifies the 
choice of the basic equations. If all the densities and velocities in equations (1) to (3) are 
divided by the reference free- stream values, these equations are derived: 
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which are the approximately correct reference expressions for boundary-layer momentum, 
displacement, and velocity thickness normalized in terms of the average or reference 
free- stream quantities. 

Often, because the differences between u^ and u§ are very small, it has been 
assumed to be justifiable without explicit substantiation, to simplify the data- reduction 
procedure or to ignore these differences in the interpretation of the data. In one proce- 
dure the boundary -layer parameters are assumed to be expressed with sufficient accu- 
racy completely in terms of local- stream- reference conditions. (See fig. 1.) This con- 
version is, of course, the correct procedure except that often the local conditions are not 
specified, and the data are presented as having been obtained under stream-flow condi- 
tions. The question thus arises as to how much of an error or difference will ensue if 
an attempt is made to correlate such data with those reported correctly. For this local- 
flow method, equations (4) to (6) convert to 
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In the other approach, it is assumed that with sufficient accuracy, p^ and Uoo 
can be substituted for p^ and u^. This procedure amounts to determining the momen- 
tum, mass-flow, and velocity losses corresponding to the Mach number losses shown by 
the shaded areas in figure 2, and neglects to isolate the changes in local velocities and 
densities induced by changes in the flow outside the boundary layer. For this free- 
stream-flow method, equations (4), (5), and (6) are modified to 
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In order to assess the possible magnitudes of the differences involved in making such 
simplifications, the following difference functions are defined: 


e a _ e A - e r _ x 

6 e* e r 

(13) 

E* - < 8 *>, - 1 
6 (5*) r 

(14) 

e A _ _ j 

(15) 

p b _ e B 1 

(16) 

e b , = < 6 *> B - 1 
6 (6*) r 

(17) 

E? = - 1 

^ (5u) r 

(18) 


where the subscript on the E defines the boundary-layer thickness parameter that is 
involved, and the superscript specifies the simplified basis on which the parameter is 
being calculated; that is, whether on local or free- stream Mach number for determining 
skin-friction momentum or other losses. Substitution of the integrals in equations (4) 
to (12) into equations (13) to (18) and the conversion of the normal distance variable y 
to y/ 5 yields 
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( 20 ) 



( 21 ) 


( 22 ) 


(23) 


(24) 


Solutions for equations (19) to (24) in closed, analytical form are obtained by 
assuming that the flow in the free stream is isentropic, and that the boundary layer has a 
power-law-type velocity profile 



(25) 


and a constant stagnation temperature equal to the free-stream value (zero heat transfer, 
Prandtl number equal to 1) which, with the aid of the thermal equation of state, yields the 
required density-velocity relationship: 
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Substitution of equations (25) and (26) into equations (19) to (24) result in (after simplifica- 
tion of the integrals, expansion, and conversion of some integrals into finite generalized 
series with a remainder term, and integration) the desired final expressions: 
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where 


a = Z-lJ:M 2 

2 6 

b = 1 + a = 1 + 2LlZ M 2 


c = Vs 
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(Equations continued on next page) 
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d = Moo 2 

e = 1 + d = 1 + M TO 2 

f = (-l) n (log e |c + l| - log e |c|) 
g = log e |c - 1| - log e |c| 
m = 1 , 2, 3, 4, * . . 

p 4{„ - i[! - 


( 33 ) 


Because of the limitations imposed by the method of evaluating some of the integrals, the 
exponential parameter n is restricted to whole numbers in equations (30) and (31). An 
illustration of the methods used in obtaining solutions is presented in appendix B. Where 
the free-stream Mach number is 0, equations (27) to (32) reduce to 
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RESULTS AND DISCUSSION 


Boundary- Layer Theory 

To assess the probable magnitudes of the differences to be expected in the use of 
the two simplified methods for reducing boundary- layer survey data to the various 
boundary- layer thicknesses, some calculations were made of the difference functions 
(eqs. (27) to (32)) for a range of the basic parameters Moo, Mg/M^, and n. The free- 
stream Mach number was varied from 0 to 10; the parameter Mg/Moo signifying the 


9 



I 

magnitude of the difference between the local and free- stream Mach numbers was varied 
from 1.00 to 1.05; and the value of the exponential factor n was varied from 6 to 14. 

An n of 6 corresponds to fully turbulent boundary layers at low Reynolds numbers and 
to values just above the conditions for transition from laminar to turbulent flow. An n 
of 14 corresponds to very high Reynolds numbers, probably on the order of 10^ or more. 

The results of the calculations where the boundary-layer thicknesses are evaluated com- 
pletely in terms of the local conditions are presented in figures 3 and 4. The results 
where the thicknesses are determined entirely in terms of the free- stream conditions 
are shown in figures 5 to 7. 

Local-flow met hod.- As expected from the formulation of the problem, the differ- 
ence functions determined by the local -flow method (eqs. (27) to (29)) are independent 
of the shape of the boundary -layer velocity profiles (figs. 3 and 4). For Mg/Moo con- 
stant and greater than 1, the difference functions are all negative at the lower free- 

A A A 

stream Mach numbers, but Eg and E g * become positive and Eg approaches 0 at 

the higher Mach numbers (fig. 3). This tendency toward more positive difference values 

with increasing Moo is generally very rapid at the lower free- stream Mach numbers 

(0.5 to 4.0), but much slower as Moo approaches 10 because of the close approach to the 

asymptotic limits existing for Moo = 00 • These trends result from the fact that Mg/M^, 

which has been taken as the basic deviation function, and the various functions of b/ e, 

which represent changes in local static pressure and local speed of sound, are opposing 

terms in isentropic flow because, physically, increases in Mg cause decreases in pg 

and the local speed of sound. As Moo approaches zero, the changes in pg and the 

local speed of sound from the assumed deviation of Mg from M m become small and 

the Mach number deviation term Mg/M^ has the dominant effect. (See eqs. (34) and 

(35) for Moo = 0.) As Moo approaches 10, the changes in pg and the local speed of 

sound resulting from the assumed deviation are very large and become the dominating 

A A 

factors in the difference terms Eg and E^*, and become essentially canceling terms 

A 

in the case of E* . 

5 u 

The results presented in figure 4 indicate that the difference terms are nearly 
linear functions of the deviation ratio Mg/M*,, within the limits of the calculations pre- 
sented herein. No calculations were made for Mg less than Moo. For equal devia- 
tions in this direction, the difference functions would be nearly identical numerically to 
those shown in figures 3 and 4, except that the positive and negative signs would be 
reversed. 

Figures 3 and 4 show that a deviation of Mg of only 1 percent (Mg/Moo = 1.01) 
will cause a difference in the calculation of the parameters 6 ^, (6*)^, and (6 U )^ of 

about -2, -1, and -1 percent, respectively, as Moo approaches 0, and will cause dif- 
ferences of about 5, 5, and nearly 0 percent, respectively, as M^ approaches 10. For 
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larger deviations the differences will be proportionately larger (within the range of cal- 
culations). For a small range of Mach number, near 1.4 for 6 ^ and 1.0 for ( 6 *)^, the 
difference in calculating these parameters is small because of the aforementioned change 
in sign of the difference functions. In view of the relatively large differences that may be 
involved over large portions of the free- stream Mach number range in reducing boundary- 
layer survey data to the various thickness parameters by this local-flow method if the 
local-flow conditions do not match the desired free-stream reference values, it appears 
desirable always to specify the local-flow conditions under which the measurements were 
made. 

Free-stream-flow method .- Where the boundary -layer thickness parameters are 
computed entirely in terms of the average free-stream conditions (free-stream-flow 
method) the values of the difference functions become dependent upon the boundary -layer 
velocity profiles (figs. 5 to 7) or on the power-profile parameter n. For M 5 /M 00 con- 
stant and greater than 1 the difference functions are again all negative at M^ near zero. 

TO TD 

As Mqq increases, E 7 ? and E* tend toward zero asymptotically, the most rapid 

v o u 

changes occurring in the Mach number range from about 0 to 4. The difference function 
B 

E * becomes positive for M^ > 1, reaches a peak value in the Mach number range from 
0 

about 3 to 4, and then declines gradually as Moo increases toward 10. 

The results presented in figures 6 and 7 indicate that the difference terms are 
nearly linear functions of the deviation ratio Mg/Moo and of the exponential parameter 
n for the range of these parameters studied herein. The differences become larger as 
n increases as a result of the relative shrinking of the low-density region near the wall. 
This decrease in the low-density region causes an increase in the mass-flow and velocity 
weighting factors involved in calculating the difference functions. 

Again, no calculations were made for Mg < Moo. For equal deviations in this 
direction the difference functions would be nearly identical numerically to those shown 
in figures 5 to 7, except that the signs would be reversed. 

An inspection of figures 5 to 7 indicates that a deviation in Mg of only 1 percent 
from the average free-stream value M^ will result in calculation differences by the use 
of the free-stream-flow technique in determining the various thickness parameters of -7 
to -16 percent at Moo = 0, -7 to 13 percent at Moo = 3, and -1 to 7 percent at M^ = 10. 
For larger deviations in Mg/M^ the differences are proportionately larger. The large 
differences that may be incurred make the use of this data- reduction method undesirable 
at any time. 

The difference curves just discussed illustrate the extreme importance of deter- 
mining and specifying Mg accurately. The assumption is made that the boundary-layer 
Mach number profile has been well defined experimentally except that an insufficient 
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number of points were taken external to the boundary layer, perhaps because of strongly 
rotational external flow, to determine Mg with great accuracy. As a consequence, the 
local Mach number external to the boundary layer was erroneously equated to Moo 
although in reality it was 1 percent higher. Then the experimenter will incur differences 
of the magnitude indicated in figures 5 to 7 although he may actually be trying to use the 
local-flow method in reducing the data. The seriousness of the problem is emphasized 
if it is recognized that even the determination of M5 to within ±1/4 percent accuracy, 
which requires sophisticated instrumentation and testing techniques, can still result in 
maximum differences up to 3 and 4 percent for some conditions. 

Experimental Wake Results 

The results of the experimental investigation are presented in table I. The basic 
data were derived from the wake surveys of reference 1. These data apply only to sta- 
tion 2 of wings F (a flat wing) and 1 (a linearly twisted wing) and only for an angle of 
attack of 0°. Station 2 was the only one at which wake surveys were made at several 
locations downstream of the wing, and restrictions of the analysis to the aforementioned 
wings at zero incidence eliminate the involvement of separated flows. Turbulent 
boundary- layer flow was assured by fixing transition near the wing leading edges. Com- 
plete details of instrumentation and testing techniques can be found in reference 1. 

Raw wake- survey data were reduced to thickness parameters by the three methods 
discussed under "Theory": the reference method, the local-flow method, and the free- 
stream-flow method. All methods, as applied, allowed for a variable M§ through the 
wake and for measured changes in local static pressure. For convenience, the reference 
method expressions for boundary-layer (and wake) momentum, displacement, and velocity 
thickness are presented in appendix C in a form most suitable for usual data- reduction 
procedures. Variable Mg values were determined by the fairing of a fictitious, linear 
pitot pressure tangential to the experimental values determined just outside either side of 
the wake. More details of this technique can be found in reference 1. Calculations indi- 
cate, however, that for the data incorporated in this report, both the effects of skewness 
of Mg and of variable local static pressure in the direction normal to the wake were 
negligible. In table I, only the momentum thicknesses have been averaged, inasmuch as 
within the simplifications incorporated into the basic reference equations (eqs. (4) to (6)), 
this parameter can remain constant as the wake expands downstream of the wing. Dis- 
placement and velocity thicknesses, however, still must decrease so that the parameters 
6 */ 0 and Sy/g approach 1 at very large distances downstream of the wing. 

Examination of the momentum thicknesses obtained by use of the reference method 
shows that the maximum scatter for the stations downstream of any one wing at one test 
Mach number was ±3.4 percent and that the average maximum scatter for all test 



conditions was on the order of ±2 percent. An analysis of all pertinent factors suggests 
that the biggest contributor to this scatter is probably the inability to determine 
with any greater accuracy. In fact, approximate calculations indicate that errors in M$ 
on the order of ±0.002 could account for nearly all the experimental differences, and it is 
doubtful that the accuracy of M5 exceeds this figure. The average maximum scatter at 
Moo = 2.01 is somewhat less than that at Moo = 1.61 and follows the trend indicated in 
figure 5, but the indication may be fortuitous because of the insufficient size of the experi- 
mental sample. The consistency of the nearly 3 -percent increase in 6 r of wing 1 over 
wing F at both test Mach numbers suggests that the indication of at least the trend may 
be meaningful, but again the sample is small. Finally, both theoretical calculations and 
the tendency of the experimental data to follow the trends derived on the isentropic-flow 
basis appear to indicate that the wing trailing- edge shock did not have a major influence 
in these tests if the wake data were properly referenced to stream-flow conditions. 

Inspection of the momentum-thickness results obtained by the local-flow method 
indicates a maximum scatter distribution nearly identical with that obtained by the refer- 
ence method. The existence of this pattern is ascribed primarily to small errors in the 
determination of M5. At M^ = 1.61, M§ is usually less than M^, and 0 A , according 
to the theoretical isentropic-flow calculations, should be slightly less than 0 r (see fig. 3 
and reverse sign), but experimentally the reverse is true. This difference in trend is 
attributed to the wing trailing-edge shock losses existing in the experimental flow and the 
relatively small isentropic-flow effects to be expected at this test Mach number. Further- 
more, the difference is so small that it is well within experimental error. At Moo = 2.01, 
6^ increases over 6 r , as expected from the theoretical isentropic-flow calculations, and 
again this difference indicates the relative smallness of nonisentropic-shock effects. 
Although the apparent difference in the use of the local-flow method in reducing the data 
of these experimental results, as compared with the reference method, is practically 
nonexistent at Moo = 1-61 and relatively small (a little over 2 percent on the average) at 
Moo = 2.01, it can be very significant at higher Mach numbers. (See fig. 3 again.) These 
differences between methods again indicate the need for specifying the flow conditions 
under which the wake data were obtained even if the difference from the stream or aver- 
age reference condition is small. 

Calculations of the momentum thickness by the free- stream-flow method result in 
very large scatter. The distribution of momentum thickness with downstream measuring 
station is no longer primarily random in character as it was for 6 r and # A , but is 
dependent upon the relative magnitude of M5 compared with Moo. The magnitude of 
the scatter is considerably larger than that expected from the boundary-layer calcula- 
tions of figure 5. This trend is due to the fact that the wake is two to three times as wide 
as a boundary layer having an equal momentum deficiency; and it can be readily inferred 
from figure 2 that as the Mach number deficiency becomes more shallow and spreads to 
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several times the boundary-layer thickness, the effects of deviations of Mg from Moo 
become proportionately more significant. Analysis also shows that as the Mach number 
deficiency shallows in the central part of the wake relative to the Mach number deficiency 
in the boundary layer close to the wall, the viscous heating or compressibility effects in 
the wake decrease relative to those computed for the boundary layer and, in effect, the 
wake difference functions become more equivalent to the boundary -layer difference func- 
tions which were computed for lower free-stream Mach numbers. For 0® this effect 
will result in increasing errors for the wake relative to the boundary layer at the same 
Moo- (See fig. 5.) The scatter in is much larger at Moo = 2.01 than at Moo = 1.61. 
This trend is primarily the result of the relatively larger deviations experienced in Mg 
at Moo = 2.01, but part of the increased scatter can also be traced to the aforementioned 
decreased compressibility effects for the wake relative to the boundary layer. This 
effect is stronger at Moo = 2.01 than at Moo = 1.61. 

At this point it is of interest to compare for this single wing test station the experi- 
mental section skin-friction coefficients derived from the average momentum thickness 
which were computed by the reference method with the theoretical values. (See table II.) 
Experimental coefficients were determined from 

Cf = 2 (^average (38) 

c s 

The theoretical values were computed by the Sommer and Short reference temperature 
method (ref. 3), and the comparison shows that the experimental skin-friction coefficients 
exceed the theoretical values by about 4 to 9 percent. A preliminary analysis indicates 
that this trend probably cannot be ascribed to effects of the transition-fixing roughness 
strip alone. 

An analysis of the displacement- and velocity-thickness results obtained in this 
experimental phase of the investigation is not presented in detail. It should suffice to 
state that this set of results conforms to the theoretical isentropic-flow boundary -layer 
calculations, with reservations for differences in wakes and boundary layers, in the same 
manner as the momentum thickness data. The only notable effect was that the decreased 
compressibility effects for the wake relative to the boundary layer resulted in consider- 
ably smaller scatter for (6*)® than might have been inferred from the theoretical 
boundary- layer calculations which neglect this effect. 

CONCLUDING REMARKS 

An analytical and experimental investigation has been made of boundary- layer and 
wake- survey data- reduction techniques in compressible flows where the local Mach 
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number just outside the boundary layer or wake does not match the average free- stream 
or normalizing value. The analytical investigation was limited to turbulent boundary 
layers with no heat transfer and to a Mach number range from 0 to 10. The experimental 
portion of the investigation was confined to zero-heat-transfer turbulent-wake surveys 
made behind a swept wing at free -stream Mach numbers of 1.61 and 2.01. 

Results of the analytical boundary-layer calculations indicate that if simplified 
data-reduction techniques that neglect the deviation of the local Mach number from the 
average free- stream or normalizing value are utilized, sizable differences in momentum, 
displacement, and velocity thicknesses are incurred over much of the Mach number range, 
for relatively small deviations in the local-flow characteristics or parameters. In gen- 
eral, use of the free-stream Mach number as a base both for determining velocity defects 
due to skin-friction and for normalizing the thickness parameters, introduces larger dif- 
ferences than use of local Mach numbers as the base for both items and interpretation of 
the results as being obtained under stream conditions. Experimental wake data indicate 
that, because of increased wake widths and reduced viscous heating or compressibility 
effects compared with boundary layers of equal momentum deficiency, the use of the free- 
stream-flow technique generally results in proportionately larger differences in wake- 
data reduction. With proper reservations for basic differences between wakes and bound- 
ary layers, the trends in the differences of the experimental wake data appear to be in 
good agreement with the analytical calculations for boundary-layer results which were 
reduced by the same method. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., October 19, 1967, 

126-13-02-11-23. 
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APPENDIX A 


CHOICE OF BASIC EQUATIONS 

The assumption that the boundary- layer momentum, mass flow, and velocity defects 
in a flow with pressure gradients can be expressed by equations (1) to (3) neglects the 
fact that the reference quantities p g and u g vary along the length of the boundary layer 
ahead of the measuring station. For thin wings the variation in the reference quantities 
is small. Furthermore, at supersonic speeds the decreases in ug and increases in pg 
experienced on the forward portions of the wing are partially compensated for by the 
reversed effects encountered on the rear parts of the wing. Thus, the overall effects are 
generally assumed to be relatively small. However, the objective of this paper is not to 
investigate this problem, but to determine the possible magnitude of differences that are 
incurred in reducing or interpreting the data, if, the boundary layer having been formed 
and measured, the reference conditions at the measuring station are misinterpreted or 
misapplied because the differences between local and free- stream conditions are believed 
to be negligibly small. For such an investigation, equations (1) to (3) should be adequate. 
For the estimation of skin friction on wings from wake surveys in pressure gradients, 
the variation of p g and u g along the wake could affect the results. Comparison of 
theoretical calculations made herein with suitable experimental results can show whether 
the omission of detailed variation of the reference quantities is of great significance. 
Simplified theoretical computations indicate that for the experimental data reported 
herein, the maximum effect is less than 2 percent and the average effect is considerably 
less. 
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APPENDIX B 


EXAMPLE OF SOLUTION OF EQUATIONS 


The equation for the difference function of the boundary- layer displacement thick- 
ness computed on the assumption that the free- stream values of p TO and u,*, can be 


substituted for Pg and Ug is 


E ?» = 


JI 


(Bl) 


/ p 6 u 5 

- pu w yN ! 

\Poo U oo 

Poo u °°/ \6/ 


With the conversion of the terms in equation (Bl) to a common denominator, the addition 
and subtraction of the term pgUg/p^u^ in the integrand of the numerator of the equa- 
tion, and the multiplication of the denominator only by the identity term PgUg/PgUg, 
equation (Bl) is transformed into 


e b * = 
6* 


I’ 

J 0 


VPco^ - P 6 U 5 + P« u 


'6"5 


pu 


p u 

CO O 


“ftH 0 


i/p 6 u g 


pu' 


Pco U o 


d(i 
0 


I, 


VP 6 u 6 V P 6 U 6 _ PU> 


Poo u c 


P 5 U 6 


d (t 


(B2) 


The integrals can now be broken up to yield 


B _ / P 0O U 0O ) 

) 


£ 4) 




(B3) 


The term p^u^/PgUg is converted into a more convenient form by the use of the 
usual (ref. 4) isentropic-flow relationships (an assumption that precludes free-stream 
shock losses) and the result is 

y+1 


Poo u oo _ 

p S u S "M 0 W 


(B4) 


where 


b=1 + ^M 2 6 
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Obviously, 


In order to evaluate the integral 



(B5) 


(B6) 


however, some assumptions have to be made about the boundary layer. It was first 
assumed that the boundary layer had a power-type velocity profile that could be repre- 
sented by a new variable u*, where 


u 


* 


u 

u 6 



(B7) 


Also, the static pressure through the boundary layer was assumed to be constant and 
equal to the value just outside the boundary layer (compatible only with zero or small 
pressure gradients) and the stagnation temperature within the boundary layer was assumed 
to be constant and equal to the free- stream value (equivalent to zero heat transfer, 

Prandtl number equal to 1). These assumptions allow the derivation, with the aid of the 
thermal equation of state, of the required density-velocity relationship 




m; 


1 - 



Substitution of equations (B7) and (B8) into equation (B6) results in 


(B8) 




du* 


(B9) 


where 


a = 



M 


2 

6 


b = l+ a = l + - Mo 

2 6 

c = 

Expansion of the integrand on the right-hand side of equation (B9) and conversion of the 
result into a finite generalized series with a remainder term that is compatible for both 
even and odd values of n result in the equation 
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c 2(m-l) u *n-2m 



(BIO) 


where 

m = 1, 2, 3, 4, . . . 



and which can now be readily integrated to yield 


where 


i 


J™_d(T) = 2| 
0 P 6 U 6 W 


>n-l 


m=p 


( f - g) - £ 


m=l 


,2(m-l) 

- 2m + 1 1 


f = (-l) n (log e |c + l| - log e |c |) 
g = log e |c - l| - log e |c| 


(Bll) 


Substitution of equations (B5) and (Bll) into equation (B3) results in the desired final 
expression 




y+1 

2(^1) 


n r 






,n-l 


C 2(m-D ' 

'^ f " Z n - 2m + 1 

m=l 


•> 


(B12) 


Because of the limitations imposed by the method of evaluating the integral of equa- 
tion (B9), the power-law parameter n is restricted to whole numbers in equation (B12). 
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APPENDIX C 


COMPUTATIONAL FORM OF THE REFERENCE BOUNDARY- LAYER- OR 
WAKE-THICKNESS PARAMETERS 

The reference expressions (where local-flow conditions at the measuring station 
do not match the desired free-stream normalizing values) for boundary- layer (and wake) 
momentum, displacement, and velocity thicknesses for conditions of no heat transfer and 
constant stagnation temperature through the boundary layer equal to the free-stream 
value are, in the form most suitable for usual data- reduction procedures, 



In equations (Cl) to (C3) p g and M g are not necessarily the local static pressure and 
local Mach number at the outer edge of the boundary layer but can be variable quantities 
if gradients in these terms are in the "potential flow" outside the boundary layer and nor- 
mal to the surface. In such instances M g normally has to be estimated by some extrap- 
olation or interpolation process as described in this paper and in reference 1. Also, the 
assumption has been made that the local static pressure has been measured and can be 
variable. For constant static pressure through the boundary layer, p g becomes con- 
stant and equal to the value ascribed to the edge of the boundary layer; for constant or 
uniform local free-stream flow, M5 becomes constant and equal to the local Mach num- 
ber just outside the boundary layer. Under these conditions the ratios involving only 
terms in p g and p M and M g and M^ can be taken outside the integral sign if 
desired. As indicated by the wake experimental results presented in the body of this 
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paper and by approximate theoretical estimates, the error in using these simplified equa- 
tions in place of the more complex and usually unsolvable fully correct equations is 
negligible over the usual range of test conditions to which they may be applied. 
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TABLE I.- CALCULATION OF WING WAKE PARAMETERS 


(a) Dimensional quantities in inches 



0.25 

.50 

.62 

.87 

i 

1.610 

1.610 

1.610 

1.610 

1.586 

1.601 

1.609 

1.591 

-1.5 

-.6 

-.1 

-1.2 

0.0398 

.0372 

.0377 

.0390 

b .0384 

b ±3.4 

0.0926 

.0828 

.0830 

.0841 

0.0545 

.0466 

.0459 

.0451 

0.0400 

.0374 

.0378 

.0391 

b .0386 

b ±3.3 

0.0920 

.0830 

.0831 

.0837 

0.0550 

.0468 

.0458 

.0456 

0.0442 

.0395 

.0377 

.0490 

b .0426 

b ±8.0 

0.0894 

.0843 

.0856 

.0744 

0.0597 

.0492 

.0459 

.0558 

Wing 1 

0.25 

1.610 

1.591 

-1.2 

0.0407 


0.0950 

0.0560 

0.0407 


0.0943 

0.0563 

0.0476 


0.0876 

0.0635 

.62 

1.610 

1.610 

0.0 

.0387 


.0855 

.0476 

.0388 


.0860 

.0476 

.0383 


.0908 

.0472 

.87 

1.610 

1.610 

0.0 

.0388 


.0841 

.0456 

.0389 


.0844 

.0456 

.0388 


.0894 

.0456 





b. 0394 

b±2.5 



b . 0395 

b ±2.4 

_i 


b .416 

b ±10.8 



Wing F 

0.25 

2.010 

2.018 

0.4 

0.0363 


0.1082 

0.0551 

0.0371 


0.1107 

0.0551 

0.0342 


0.1337 

0.0528 

.62 

2.010 

2.036 

1.3 

.0359 


.1020 

.0475 

.0367 


.1050 

.0471 

.0289 


.1351 

.0395 

.87 

2.010 

2.055 

2.2 

.0364 


.1022 

.0460 

.0374 


.1063 

.0455 

.0242 


.1368 

.0342 





b. 0362 

b±0.6 



b . 0371 

b ±0.9 



b.0291 

b ±26.0 



Wing 1 

0.25 

2.010 

2.089 

0.9 

0.0379 


0.1125 

0.0560 

0.0383 


0.1142 

0.0558 

0.0313 


0.1331 

0.0488 

.62 

2.010 

2.046 

1.8 

.0363 


.1039 

.0484 

.0370 


.1067 

.0479 

.0287 


.1284 

.0394 

.87 

2.010 

2.064 

2.7 

.0374 


.1053 

.0471 

.0381 


.1089 

.0465 

.0156 


.1519 

.0253 





b.0372 

b ±2.1 



b .378 

b±2.0 



b .252 

b ±33.5 




co 

00 


^Scatter is based on mean of high and low values of 6. 
^Average values. 
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TABLE I.- CALCULATION OF WING WAKE PARAMETERS - Concluded 


(b) Dimensional quantities in centimeters 



u — 

Wake parameters 


1 

VActCII 11 

uiiiutri 

— 
Reference method 

Local-flow method 

Free-stream method 

Station, 

x/c s 

Moo 

Mg 

Difference, 

percent 

e r , 

cm 

Scatter 

in 6 V , 
percent 

(a) 

(6*) r , 

cm 

( 5 u) r > 

cm 

e A , 

cm 

Scatter 

in 0 A , 
percent 

(a) 

(6*) A , 

cm 

; 

( 6 u) A > 

cm 

e B , 

cm 

Scatter 

in 0 B , 
percent 

(a) 

(5*) B , 

cm 

(6u) B , 

cm 


Wing F 


0.25 

1.610 

1.586 

-1.5 

0.1011 


0.2352 

0.1384 

0.1016 


0.2337 

0.1397 

0.1123 


0.2271 

0.1516 

.50 

1.610 

1.601 

-.6 

.0945 


.2103 

.1184 

.0950 


.2108 

.1189 

.1003 


.2141 

.1250 

.62 

1.610 

1.609 

-.1 

.0958 


.2108 

.1166 

.0960 


.2111 

.1163 

.0958 


.2174 

.1166 

.87 

1.610 

1.591 

-1.2 

1 .0991 

.2136 

.1146 

.0993 

i 

.2126 

.1158 

.1245 


.1890 

.1417 


b .0975 b ±3.4 b .0980 b ±3.3 b .1082 b ±8.0 


Wing 1 


0.25 

1.610 

1.591 

-1.2 

0.1034 


0.2413 

0.1422 

0.1034 


0.2395 

0.1430 

0.1209 


0.2225 

0.1613 

.62 

1.610 

1.610 

0.0 

.0983 


.2172 

.1209 

.0986 


.2184 

.1209 

.0973 


.2306 

.1199 

.87 

1.610 

1.610 

0.0 

.0986 


.2136 

.1158 

.0988 


.2144 

.1158 

.0986 


.2271 

.1158 





b .1001 

b ±2.5 



b .1003 

' b ±2.4 


i 

i 

b .1056 

b ±10.8 




Wing F 


0.25 

2.010 

2.018 

0.4 

0.0922 


0.2748 

0.1400 

0.0942 


0.2812 

0.1400 

0.0869 


0.3396 

0.1341 

.62 

2.010 

2.036 

1.3 

.0912 


.2591 

.1207 

.0932 


.2667 

.1196 

.0734 


.3432 

.1003 

.87 

2.010 

2.055 

2.2 

.0925 


.2596 

.1168 

.0950 


.2700 

.1156 

.0615 


.3475 

.0869 


i 



b .0919 

b ±0.6 



b .0941 

b ±0.9 



b .0739 

b ±26.0 




Wing 1 


0.25 

2.010 

2.029 

0.9 ! 

i 

0.0963 

1 

0.2858 ' 

0 . 1422 ! 

0.0973 ' 


0.2901 

0.1417 

0.0795 


0.3381 

0.1240 

.62 

2.010 

2.046 

1.8 

.0922 


.2639 

.1229 

.0940 


.2710 

.1217 

.0729 


.3261 

.1001 

.87 

2.010 

2.064 

2.7 

.0950 


.2675 

.1196 

.0968 


.2766 

.1181 

.0396 


.3858 

.0643 





b .0945 

b ±2.1 



. b .0960 

b ±2.0 



b .0640 

b ±33.5 




a Scatter is based on mean of high and low values of 6 . 
b Average values. 



TABLE II.- COMPARISON OF CORRECT EXPERIMENTAL 
SKIN FRICTION WITH THEORY 


Wing 

Mach 

number 

Section skin-friction coefficient, 

Cf 

Difference, 

percent 


Experimental 

Theoretical 

F 

1.61 

0.00610 

0.00587 

3.9 

1 

1.61 

.00625 

.00587 

6.5 

F 

2.01 

.00573 

.00540 

6.1 

1 

2.01 

.00590 

.00540 

9.3 
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